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Abstract 



A review of methods for finding general expressions for matrix el- 
ements (non-diagonal with respect to configurations included) of any 
one- and two-particle operator for an arbitrary number of shells in an 
atomic configuration is given. These methods are compared in various 
aspects, and the advantages or shortcomings of each particular method 
are discussed. Efficient method to find the abovementioned quantities 
in LS coupling is presented, based on the use of symmetry properties 
of operators and matrix elements in three spaces (orbital, spin and 
quasispin), second quantization in coupled tensorial form, graphical 
technique and Wick's theorem. This allows to efficiently account for 
correlation effects practically for any atom and ion of periodical table. 
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1 Introduction 



Modern atomic spectroscopy studies the structure and properties of prac- 
tically any atom of the periodic table as well as of ions of any ionization 
degree. Particular attention is paid to their energy spectra. For the in- 
vestigations of many-electron atoms and ions, it is of great importance to 
combine experimental and theoretical methods. Nowadays the possibilities 
of theoretical spectroscopy are much enlarged thanks to the wide use of 
powerful computers. Theoretical methods utilized must be fairly universal 
and must ensure reasonably accurate values of physical quantities studied. 

Many-electron atom usually is considered as many-body problem and is 
described by the wave function constructed from the wave functions of one 
electron, moving in the central nuclear charge field and in the screening field 
of the remaining electrons. Then the wave function of this electron may be 
represented as a product of radial and spin-angular parts. The radial part is 
usually found by solving various modifications of the Hartree-Fock equations 
and can be represented in a numerical or analytical forms (Froese Fischer pQ) 
whereas the angular part is expressed in terms of spherical functions. Then 
the wave function of the whole atom can be constructed in some standard 
way (Cowan 2 , Jucys and Savukynas [3], Nikitin and Rudzikas fU) start- 
ing with these one-electron functions and may be used further on for the 
calculations of any matrix elements representing physical quantities. 

For obtaining the values of atomic quantities it is necessary to solve 
so-called eigenvalue problem 



where is the wave function of the system under investigation and H is 
its Hamiltonian. In order to obtain accurate values of atomic quantities 
it is necessary to account for relativistic and correlation effects. It turned 
out that for very large variety of atoms and their ionization degrees the 
relativistic effects may be taken into account fairly accurately as Breit-Pauli 
corrections (Nikitin and Rudzikas [I], Rudzikas 5 ). It is convenient to 
present the Hamiltonian as consisting of two parts in Breit-Pauli approach, 
namely, 



where Hnr is the non-relativistic Hamiltonian and Hr stands for the rela- 
tivistic contibution. Non-relativistic Hamiltonian is the sum of the kinetic 




Hbp = Hnr + Hr 



(2) 
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energy of electrons T, their potential energy P, and the electrostatic electron 
interaction Q 

H NR = T + P + Q. (3) 

The first two terms, T and P are the one-particle operators, whereas the 
third term, Q, is two-particle operator. The Hr may also be subdivided 
into non-fine structure and fine structure contributions 

Hr = H Non-Fine + Hpine- (4) 

The non-fine structure contributions 

HNon-Fine = Huass + Hoi + Hd2 + HssC + Hoo (5) 

shift non-relativistic energy levels without any splitting of them. The mass- 
velocity term H^ass describes the variation of mass with velocity. The one- 
and two-body Darwin (contact) terms Hm and B.r>i are the corrections of 
the one-electron Dirac equation due to the retardation of the electromagnetic 
field produced by an electron (contact interaction). The spin-spin-contact 
term H$sc accounts for the interaction of the spin magnetic moments of 
two electrons occupying the same space. The orbit-orbit interaction Hoo 
accounts for the interaction of two orbital moments. 
The fine-structure contributions 

HFine = Hso + HsoO + Hss (6) 

split the non-relativistic energy levels (terms) into a series of closely-spaced 
levels. The most important of these is the nuclear spin-orbit interaction 
Hgo (spin-own-orbit) representing the interaction of the spin and angular 
magnetic moments of an electron in the field of the nucleus. The spin-other- 
orbit Hsoo an d spin-spin Hss contributions in a rough sense may be viewed 
as corrections to the nuclear spin-orbit interaction due to the presence of the 
other electrons in the system. 

Unfortunately, practical calculations show that all realistic atomic Hamil- 
tonians do not lead straightforwardly to eigenvalue problem (^Q) . Actually we 
have to calculate all non-zero matrix elements of the Hamiltonian considered 
including those non-diagonal with respect to electronic configurations, then 
to form energy matrix, to diagonalize it, obtaining in this way the values 
of the energy levels as well as the eigenfunctions (the wave functions in the 
intermediate coupling scheme). The latter may be used then to calculate 
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electronic transitions as well as the other properties and processes. Such a 
necessity raises special requirements for the theory. 

The total matrix element of each term of the energy operator in the case 
of complex electronic configuration will consist of matrix elements, describ- 
ing the interaction inside each shell of equivalent electrons as well as between 
these shells. Going beyond the single-configuration approximation we have 
to be able to take into account in the same way non-diagonal, with respect 
to configurations, matrix elements. 

To find the expressions for the matrix elements of all terms of the Hamil- 
tonian considered for complex electronic configurations, having several open 
shells, is a task very far from the trivial one. A considerable part of the ef- 
fort must be devoted to coping with integrations over spin-angular variables, 
occurring in the matrix elements of the operators under consideration. This 
paper presents the general methodology, leading to optimal expressions for 
operators and matrix elements. 

A number of methodologies to calculate the angular parts of matrix 
elements exists in literature. Many existing codes for integrating the spin- 
angular parts of matrix elements (Glass (Hj, Glass and Hibbert [Jj, Grant [H], 
Burke et al (§]) are based on the computational scheme proposed by Fano |1U| 
This methodology is based on having the total wavefunction of an atom built 
from the antisymmetrized wavefunctions of separate shells, and this antisym- 
metrization is done via coefficients of fractional parentage. The shells are 
coupled one to another via their angular momenta. So, the finding of matrix 
elements amounts to finding the recoupling matrices and the coefficients of 
fractional parentage. 

Suppose that we have a bra function with u shells in LS coupling: 

^™ (LSM L M S ) = 
= (n 1 lin 2 l 2 ...n u l u a 1 L 1 Sia 2 L 2 S 2 ...a u L u S u ALSM L M s \ (7) 

and a ket function: 

ifi* (L'S'M' L M' S ) = 
= \n 1 l 1 n 2 l 2 ...n u l u a' 1 L' 1 S' 1 a' 2 L' 2 S 2 ...a' u L' u S' u A' V S' M' L M' S ), (8) 

where A stands for all intermediate quantum numbers, depending on the 
order of coupling of momenta LiS{. Using the Wigner-Eckart theorem in 
LS space we shift from the matrix element of any two-particle operator G 
between functions ® and (jHJ) to the submatrix element (reduced matrix) 

(V£ ra (LS) \\G( K1K2k ' a ^\\ip^ et (L'S')) of this operator. When the two- 
particle operator acts upon four distinct shells, then the finding general 
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expressions of matrix elements, according to methodology by Fano ^U], is 
based upon the formula 



(lpl ra (LS) ||($(«i*2Avxi<X2fc)||^et (US')) ~ 

£ (_i)A [jy. (iV . _ ^ (jy/, _ ,5(/,i'))] 1/2 x 

7ij Aj ,rij Aj ,71^ A ; / ,riji\j/ 

x (zf* c^SiHZf*- 1 (a^),fc) (Zf aj-L^HZf*" 1 (n^V) ./ ; ) x 



[R d (\ U Xj, X'i, Xj, Ki, K2, k, (71,0-2, o-, A bra ,A kl 
x (1 + tffrjW.i'))" 1 (n.A^A.H^^^^^IKA^A;. 



-i? e ( Xi , Xj , X[ , X'j, Ki , K 2 , K, O-i , <7 2 , (7, A bra , A ket 

x (l + 5(i,j))(l -<5(i',/)) (n 4 A f n i A i ||^' t »*.^*)||Ti^A;.n^)} 

(9) 

where A is a phase factor, (see for example in jS]), 

A bra = (LiSi, LjSjiL'iS'i, L'jS'^j" a is the array for the bra function shells' 

terms, and similarly for A ket . The coefficient (l N aLS\\l N ~ l {a'L'S')l\ is 
a fractional parentage coefficient, and coefficients 
R d (x i ,X j ,X' i ,X' j ,K 1 ,K2,K,a 1 ,a 2 ,cr,A bra ,A ket ^ and 

\Xi, Xj, X^, Xj,Ki,K2, K,ai,a2, a, A bra , A ket ^j are the recoupling matri- 
ces in I and s - spaces of direct and exchange terms, respectively. For more 
detailes on recoupling matrices see in Grant [Sj, Burke el al [SJ. 

The summation in expression Q implies that the quantum numbers n, 
X of all participating shells are included. There are four such pairs of n, X 
in the sum. 

In essence, the Fano calculation scheme consists of evaluating recou- 
pling matrices. Although such an approach uses classical Racah algebra 

1121 I14j on the level of coefficients of fractional parentage, it may be 
necessary to carry out multiple summations over intermediate terms. Due 
to these summations and the complexity of the recoupling matrix itself, 
the associated computer codes become rather time consuming. Jucys and 
Vizbaraite |15j proposed to use the two-electron coefficients of fractioanl 
parentage instead of ordinary ones, in matrix elements' calculations, but 
even that did not solve the abovementioned problems. A solution to this 
problem was found by Burke et al [S]. They tabulated separate standard 
parts of recoupling matrices along with coefficients of fractional parentage 
at the beginning of a calculation and used them further on to calculate the 
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needed coefficients. Computer codes by Glass [H], Glass and Hibbert [7], 
Grant jS], Burke et al utilize the program NJSYM (Burke (THj) or NJ- 
GRAF (Bar-Shalom and Klapisch |17| ) for the calculation of recoupling ma- 
trices. Both are rather time consuming when calculating matrix elements of 
complex operators or electronic configurations with many open shells. 

In order to simplify the calculations, Cowan [2] suggested that matrix 
elements be grouped into "Classes" (see Cowan (2] Figure 13-5). Unfor- 
tunately, this approach was not generalized to all two-electron operators. 
Perhaps for this reason Cowan's approach is not very popular although the 
program itself, based on this approach, is widely used. 

Many approaches for the calculation of spin-angular coefficients (Glass (Hj, 
Glass and Hibbert ,7j, Grant |H|, Burke et al [S]) are based on the usage of 
Racah algebra only on the level of coefficients of fractional parentage. A few 
authors (Jucys and Savukynas 3 , Cowan |2|) utilize the unit tensors, sim- 
plifying the calculations in this way, because use can be made of the tables 
of unit tensors and selection rules can be used prior to computation to check 
whether the spin-angular coefficients are zero or not. Moreover, the recou- 
pling matrices themselves have a simpler form. Unfortunately, these ideas 
were applied only to diagonal matrix elements with respect to configurations, 
though Cowan [2j suggested the usage of unit tensors for non-diagonal ones 
as well. 

All the above mentioned approaches were applied in the coordinate repre- 
sentation. The second quantization formalism (Judd [181 119j . Rudzikas and 
Kaniauskas [201 and Rudzikas f2|) has a number of advantages compared to 
coordinate representation. First of all, it is much easier to find algebraic 
expressions for complex operators and their matrix elements, when rely- 
ing on second quantization formalism. It has contributed significantly to 
the successful development of perturbation theory (see Lindgren and Morri- 
son |25) Merkelis et al (22]) and orthogonal operators (Uylings [2S])) where 
three-particle operators already occur. Uylings |21] suggested a fairly simple 
approach for dealing with separate cases of three-particle operators. 

Moreover, in the second quantization approach the quasispin formalism 
was efficiently developed by Innes (23 , Spakauskas et al (2H1 E2 , Rudzikas 
and Kaniauskas |2Uj . Fano and Rau |28|. The main advantage of this ap- 
proach is that applying the quasispin method for calculating the matrix 
elements of any operator, we can use all advantages of the new version of 
Racah algebra (see Rudzikas p] ) for integration of spin-angular part of any 
one- and two-particle operator. For example, the reduced coefficients of frac- 
tional parentage are independent of the occupation number of the shell (see 
Rudzikas and Kaniauskas [2111 > Gaigalas et al |2H| ) • All this enabled Merkelis 
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and Gaigalas 30 to work out a general perturbation theory approach for 
complex cases of several open shells. In the paper by Merkelis [3^ a detailed 
review of a version of graphical methodology is presented that allows one to 
represent the operators graphically and to find the matrix elements of these 
operators using diagrammatic technique. 

The majority of methods and computer codes of finding angular coeffi- 
cients discussed above were faced with a number of problems, the main of 
these being: 

• The high demand of CPU time for calculating the angular parts of ma- 
trix elements even on the modern computer. Therefore the high accu- 
racy of characteristics of atomic quantities somtimes is even unattain- 
able. 

• The methods are applicable in practice only for the comparatively 
simple systems, because in the programs based on classical Racah 
algebra the treatment of recoupling matrices is rather complicated, 
especialy when finding the matrix elements of Breit-Pauli operators 
between complex configurations. 

• The configurations with open /-shells must often be included in theo- 
retical calculations. This causes problems in a number of methodolo- 
gies, because the complete account of /- shells implies using a large 
number of coefficients of fractional parentage. 

Gaigalas and Rudzikas [321, Gaigalas et al 33 j, Gaigalas et al 
Gaigalas et al [22] and Gaigalas and Rudzikas jSS] suggested an efficient 
and general approach for finding the spin-angular parts of matrix elements 
of atomic interactions, relying on the combination of the second-quantization 
approach in the coupled tensorial form, the generalized graphical technique 
and angular momentum theory in orbital, spin and quasispin spaces as well 
as on the symetry properties of the quantities considered. This approach is 
free of previous shortcomings. 

It is the main goal of this work to present all this methodology consis- 
tently and in a unifying manner, paying the special attention to its main 
ideas. Also, we aim at a detailed discussion of obtaining the efficient ten- 
sorial expressions of a two-particle operator, as well as the analytical ex- 
pressions for recoupling matrices. In addition to that, in this work we aim 
at comparing this methodology to other calculation of angular coefficients 
schemes. 
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Prom what is said above we see that the treatment of angular parts of 
matrix elements is a many-sided question. The classical scheme by Fano was 
developed in various aspects. The methods were developed that used Racah 
algebra |llU12llT3*lll4| on a higher level, or used the new version of the Racah 
algebra (angular momentum theory in orbital, spin and quasispin spaces) 
(see Rudzikas jUSl)- The ways were searched for, to obtain the expressions 
that had simpler recoupling matrices. In order to discuss all that in more 
detail, to compare the existing methods of obtaining the angular parts and 
to mark the advantages of one method or another, the finding angular parts 
must be looked at from different angles. 

Thus, the first thing to discuss is the expression for any physical op- 
erator (see sections 2,3,4), because already the form of it determines the 
level of application of tensor algebra in the calculation of matrix elements 
of this operator. Also, it is very important to make clear to what extent 
the Racah algebra is exploited in one or another methodology of angular 
parts treatment. This helps to mark the advantages of one method against 
another (see section 5). In addition, it is of importance to compare the ways 
of calculating the recoupling matrix in various methods (see section 6). 



2 Tensorial Expressions for Two-particle Opera- 
tors 

It is well known in the literature that a scalar two-particle operator may 
be presented the following tensorial form (see Jucys and Savukynas [3], 
Glass 0): 



frwiwo*) = J-g ( n , rj ) £ (-l) k - p [# 1<Tl) x t 2a2) ] { l^ , (10) 
i>j v 



p-p 



where g {ri,rj) is the radial part of operator, g{ Kiai ^ i s a tensor acting upon 
the orbital and spin variables of the i-th function, k%, K2 are the ranks of 
operator acting in orbital space, and oi, 02 are the ranks of operator acting 
in spin space. 

All the above mentioned approaches were usually applied in the coordi- 
nate representation. Now we will investigate into the second quantization 
formalism, which is broadly applied in atomic physics as well. 

A two-particle operator in second quantization method is written as 
follows: 
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G (W,«™*0 = G(ij,i'j>) = - £ a iaj a],4 (i,j\g\i',f) , 

(11) 

where i = nilism^nis^ (i,j\g\i',j') is the two-electron matrix element of 
operator Q( K i K 2k,aia2k) ^ anc j a . j g ^ e e i ec ^ ron creation and aj electron anni- 
hilation operators. Meanwhile two tensorial forms are well known in second 
quantization. In the first form the operators of second quantization follow 
in the normal order: 



giMkwk) = ^2 G I {ij,i'j') = 

nili,n<ilj,n f .l f .,n f .l f - 

Yl J2 J2 {-\) k ~ p [Ki2,/ci 2 ,ffi2, 0-12] 1/2 



x (n i Ain i A i ||c/( Kl ' t2fc ' <T1,J2fe )||n^A / i n / 7 A / 7 .) x 
h lj K \2 ] j s 

K\ K2 k 

h Ij «12 J [ s S Oyi 
(^120-12) T~fA') ~(A' 



j(Ai) x a (X 3 ) 



x a 



12°12 J 



(fcjfc) 



P-P 



(12) 



where [a, 5] = (2a + 1) (26+1), A = Is, (n i A < n J -A i ||^* ,1 ' ,afc '« T1 ^ fc )||n^n^A^ 
is the two-electron submatrix (reduced matrix) element of operator q( k i k 2k,^i^20-) 
and tensor a is defined as (see for example Rudzikas [3]) 



l m\~- 



rn x ■ 



(-1) "a 



(13) 



The product of tensors 



a (AJ x a (A J j 

l(KlK2fc,0"10"2fc) 



x a 



(fcjfc) 



P-P 



denotes tensorial part of operator G\ 

In another form the second quantization operators are coupled by pairs 
consisting of electron creation and annihilation operators. In coupled ten- 
sorial form: 
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a 



II 



E G n (ij,i>j>) 



ndi.ndj.n'-V-.n'-l'- 



E E (-!) fc " P (nA^A.Hs^^^^lK^n^^) x 



ndi-nd-j ,n .l . ,n l . V 
tt ' jj' i?' j j 



[«1,«2, 0-1,0-2 



-1/2 



.(AO 



x a 



ft) 



(fClCTl) 



("I) 



«1 
ft 




,(Ai) 



x a 



X) 



(fcfe) 



(fefe) 



5 (rijljin'jl 



n 'i) 



p-p 



(14) 

The expression ()12[) consists of only one tensorial product whereas (|14|) 
has two, but the summation in the first formula is also over intermediate 
ranks k\2, (J12, n' l2 and a' 12 , complicating in this way the calculations. The 
advantages or disadvantages of these alternative forms of arbitrary two- 
electron operator may be revealed in practical applications. 

In these forms the product of second quantization operators denotes 
tensorial part of operator Q( K i K 2k,cria2k) _ p or i ns t anC g ; the tensorial structure 
of electrostatic (Coulomb) electron interaction operator is K\ = K2 = k, o\ = 
o<i = (Jucys and Savukynas [3|), and only the two-electron submatrix 
elements (jiiXiUj Xj\\g^ K1K2k ' air72k ^\\n' i X > i n'j X'^j of these operators are different. 
In the case of electrostatic interaction: 



(m Xirij Xj I \g { couS I \ n 'iK n 'j Aj) = 
= 2[kf 2 (h\\C^\\l^ (ijWC^Wl'j) R k {rnhn&njljnft) . (15) 

Prom Q15j) , by (|12[) and (|14|) , we finally obtain the following two secondly 
quantized expressions for Coulomb operator (see Merkelis et al |37j): 



Vi 



i E E (-i)^ +fe+Kl2 [-i2^i 2 ] 1/2 (k\\cw\\ii) x 

nilinjljn'IWjl'jKizcT 12k 



(Ai) x a (Aj) 



(Kl2fl2) 



(KI2CT12)' 



(00) 



x a 



(16) 
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= E E (fcllC^llU) (hWC^Wl'j) Rk (nMlln^n'/j) x 



nilirijLn'-l'-n '-/' A; 



x W 



-1/2 



x a 



a^x a 



(00) 



+ 



+(2[yy 



-1/2 



a^x a 



( a ;) 



(00) 



(17) 

The tensorial expressions for orbit-orbit and other physical operators in 
second quantization form may be obtained in the same manner. 

It is worth mentioning that the expressions ()16|) and (|17|) embrace, al- 
ready in an operator form, both the diagonal interaction terms, relative to 
configurations, and the non-diagonal ones. Non-diagonal terms define the 
interaction between all the possible electron distributions over the config- 
urations considered, differing by quantum numbers of not more than two 
electrons. 

The merits of representing operators in one form or another (|lfi|) or (|17|) 
are mostly determined by the technique used to find their matrix elements 
and quantities in terms of which they are expressed. 

In the paper Gaigalas and Rudzikas |32| it was shown that the tensorial 
forms (flT)]). (fT2|) . (|14j) of two-particle operator do not take full advantage of 
tensor algebra. The most characteristic examples are when configurations 
considered have many open shells, or when the non-diagonal matrix elements 
are seeked. 

In the paper Gaigalas et al jHSj the following optimal expression of two- 
particle operator is proposed, which allows one to make the most of the 
advantages of Racah algebra (see Racah 1121 fTTH 114j ). 



K12,fl2Ai2' <T 12 



5(u,2) 
5(u,3) 



(«120-12) 



+E [ D( - laS) x D^i 33 ) 
+^f[l)M x nMl (KlM2) 



(kk) 

P-P 
1 (kk) 



P-P 

(kk) 
P-P 



5(u,4) 



(18) 

Whereas in traditional expressions, e. g. ()12|) . the summation runs 
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over the principle and the orbital quantum numbers of open shells without 
detailing these, in the expression written above the first term represents the 
case of a two-particle operator acting upon the same shell n a X a , the second 
term corresponds to operator G^ KlK2k ' air72k ' > acting upon two different shells 
n a X a , np\p. When operator Q( K i K -2 k ^i (T 2 k ) ac ts upon three shells the third 
term in (|18|) must be considered and when it acts upon four - the fourth 
one. We define in this expression the shells n a X a , npXp, n 7 A 7 , n$X$ to be 
different. 

The tensorial part of a two-particle operator is expressed in terms of oper- 
ators of the type (nA, S), B( kk \nX,Z), C^ kk \nX, S), D^ ls \ E^ kk \nX, 3). 
They correspond to one of the forms: 

< A) , (19) 



a m q i X a m q l 



(reici) 



(20) 



a m q \ X 



a m q \ X a mJ 



(kicti) 



(K2CT2) 



(21) 



x a 



(?A) 
m qS 



(22) 



*m 9 i 



,(<?A) 
"1^2 



(kicti) 



,(9A) 
rn q 3 



,(?A) 



(K2CF2) 



{kk) 



(23) 



For example, if we take a two-particle operator acting upon two shells, 
then we see from expression (|18|) that the angular part of two-particle oper- 
ator is expressed via operators B^ Kl2<Jl2 ^ (n a A a ,H) and C^ 12 " 12 ) (ngA^E). 
In a case when the operator Q^i^k^^k) ac ^ s - n a manner that two 
operators of second quantization act upon one shell and two act upon an- 
other, the J B( K12CT12 ) (n Q A Q ,H) and cK 2 ^ 2 ) (n p Xp,E) are expressed as CT . 
But in a case when three operators of second quantization act upon one shell 
and one acts upon another, then B^ 12(712 ^ (n a X a ,E) and C^' 12 ' 7 ") (ngA^H) 
are expressed either as (|19|) and (|21|) or (|19|) and (|22|). 

In writing down the expressions ()19|) - Q23|) the quasispin formalism was 

used, where am\ and a m ^ are components of the tensor , having in 



additional quasispin space the rank q = \ and projections m q = ±i, i.e. 



a<* A > 



""miin s 



(24) 
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and 



^mim s 



(25) 



In the expression (|18j) u is the overall number of shells acted upon by 
a given tensorial product of creation/annihilation operators. Parameter S 
implies the whole array of parameters (and sometimes an internal summation 
over some of these is implied, as well) that connect the amplitudes of 
tensorial products of creation/annihilation operators in the expression (|18|) 
to these tensorial products (see Gaigalas et al (HS1)- Also, attention must 
be paid to the fact that the ranks K\, K2, k, o~\, o~i and a are also included 
into the parameter S. The amplitudes (H) are all proportional to the 
submatrix element of a two-particle operator g, 



(H) ~ (niXiUjXj \\g\\ rtyAj/n^Aj/) 



(26) 



To obtain the expression of a concrete physical operator, analogous to 
expression ()18(l , the tensorial structure of the operator and the two-particle 
matrix elements ()26|) must be known. We shall investigate this now. 

The electrostatic ( Coulomb ) electron interaction operator jj Coulomb itself 
contains the tensorial structure 



H 



Coulomb 



(fcfcO,000) 
Coulomb 



(27) 



and its submatrix element is 



H. 



(fcfcO.OOO) 
Coulomb 



n^i Xi'Tlji Xji^J 



2[k]V 2 (k\\cM\\li>) (lj\\c^\\l f )R k ( 



(28) 



The spin-spin operator H ss itself contains tensorial structure of two dif- 
ferent types, summed over k (Gaigalas and Rudzikas |35j). 



jjss = ^2 \jj(k+lk-12,ll2) _|_ ^-(fc-ife+12,112) 



(29) 



Their submatrix elements are (Jucys and Savukynas [H]) 



( 



UjXiHjXj 



H. 



(fc+ifc- 12,112) 



ni/Xi'Uj'Xj'^j = -^\J (2k + 3) 



(5). 



x(li C {k+1 ^ h')(lj C {k -^ If)^- 1 (niUnjljM'U'nj'lj') 



(30) 



13 



Er (fc-lJfc+12,112) 



nitXitUj/XA = ^W(2/c + 3) {5) x 



(niXirijXj 

x (k ||c ( * -i) || k) (lj ||c {fc+1) || lj) N k - 1 (njljnih^flfnffli') , (31) 

where we use a shorthand notation 

(2k + 3) {5) = (2k + 3) (2k + 2) (2k + 1) (2k) (2k - 1) and radial integral 
(GOJ), (|3*T|) is defined as in Glass and Hibbert [7j: 



N k (mlinjlj,ni'li'rij>lj') 

X /o°° Jo°° Pi (n) Pi (r 2 ) 4ae(ri - r 2 )P v (n) P,v (r 2 ) dndr 2 , (32) 

' i 



where e(x) is a Heaviside step-function 

e(x) = 



1; for x > 0, 
0; for x < 0. 



(33) 



The spin- other- orbit operator fP°° itself contains tensorial structure of 
six different types, summed over k (see Gaigalas et al |34j ) : 



jjsso — \jj(k-lkl,101) , |j(*;-lfcl,011) , 
/ ^ sso sso 

k 

rr{kkl,lQX) rr(fefcl,011) rr(fe+lfcl,101) ^-(fc+lfcl.Oll) 
Usso T tLsso \ -Hsso T J^sso 



(34) 



Their submatrix elements are: 





Tr{k — lfcl,(Ti(T2l) 


n 


i' Xi'Tlji Xji i 


x (Zi + Zj/ 


- ]fe + l)(fc-J 


i + U') (k + k 


x (k)- 1 / 2 








h 


c(*) 





(niAinjAj |pJoo 1, ' T1<T2l) || Tii/Ai/n^Ajv) = -2 • 2 CT2 (2ifc + 1) 1/2 (/, \\c {k) \\ h>) 
x (/j ||c7( fc )|| l r ) [(k(k + 1))~ 1/2 (k (k + 1) - k (k + 1) - (Z;/ + 1)) 
x {(& + 1) N k ~ 2 (njljUilijUj/ljini'lii) — kN k (nilinjlj,ni>lirnj>ljt)\ 
-2 (k (k + 1)) 1/2 V k ~ l (nilirijlj^idi'njdjr)} , 

(36) 
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(niXiUjXj ||# s ( o^ lfcl ' CTl<T2l) || ni'Xi'UfXj^ = 2 • 2< 72 {(2k + 1) (2k + 3) 
x (k + l v -k)(k-U + l v + l)(k + h- h> + l)(k + h + l v + 2)} 1/2 
x (k + iy 1/2 (li \c [k) \ k>) (lj \\c {k) I lA N k (mknjlj^idimjdj,) . 

The radial integrals in (|33|) - (|3*7j) are (see Glass and Hibbert [7]): 



(37) 



V k (nilinjlj,ni'liinjilj>) 
£ /o°° Jo°° ^ (n) ^ (r 2 ) Cr^P. (n) (r 2 ) d n dr 2 . 



(38) 



The tensorial form of orbit-orbit operator is (see Eissner e£ al |38j ^ 



P 



(fcfcO,000) 
ool 



+p, 



(fcfcO,000) 
oo2 



+ P 



(fcfcO.OOO) 
oo3 



+ P 



(fcfcO,00)\ 
oo4 J 



(39) 



P. 



The sum of submatrix elements of three terms (^fco.ooo) ^ ^(fefco,ooo) an( ^ 

(fcfcO,000) . 



oo 1 



is equal to (see Badnell 39 ): 



Tl^XiTlj Xj 



P 



(A;fcO,000) 
ool 



+ P 



(fcfcO.OOO) 
oo2 



+ P 



(A:fcO,000) 



oo4 



flj' Xi'Tlji Xj/^j 



= j — 2A; (/c + 1) ^T fe+1 (nilinjlj^ni'liirijilji) — T k 1 (riilirijlj, nidi>nj/lji)^j 
+ (ii (|< + 1) - A; (* + 1) - ^ & + 1)) 
x (u k+1 (mlinjlj,ni/liinj/lji) — U k ~ l (riilirijlj, ni'li/rijdj/) 

+ (lj (lj + 1) - k (k + 1) - lj, (lj, + I)) 
x (ll k+1 (rijljnili,njiljini/lii) — U k ~ 1 (njljnili,njdj>ni'li>) 
+\ (k (k + 1) - k (k + 1) - h> (k> + 1)) (l 3 - (lj + 1) - fc (k + 1) - j/^/ + 1)) 
x fc(2fc-i) (^ r/c_2 (nikrijlj,ni'li'rij'lj>) + iV fc ~ 2 (rijljnili,njdj>nidi>)) 
- ( fc+ i fc )(2fc+3) (niknjlj,ni>li>nj'lji) + iV fe (njljml^njdjinidii)^ } 

x(l-<S(fc,0)), 

(40) 

where 

T fc (nilinjlj,ni>li/njilj>) = 4( - 2 ^. +1 - ) 

x Jo 00 /o°° P (»l) P, (-2) (& + i) Pi' (n) (£ + i) (r 2 ) dr lf *r 2 , 

(41) 
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U (nilinjlj,rii'li'nj'lj') = ^k+i) 



x / °° Jo°° Pi (n)Pj fa) (* " l)4fec(n - r 2 ) -(k + 2) Wr 2 - n) 

\ r l r 2 

*Pi> (n)(4 + 4)Pi< (r 2 )d ri dr 2 . 

(42) 

The submatrix element of remaining term H^ ' is: 

(niXiUjXj #^0,000) m/Xi'TifXA = 2^/WTT^ ! ^ T j 
x + ij/ + k + 2) (Zi + l v - k) (k -l v + k + l) (l v -h + k + 1) 
x (lj + l f + k + 2) x (Zj + lj> - k) (lj - If + k + l) (If -lj + k + l)) 1/2 
x(z i ||C( fc + 1 )||^)(/ i ||C( fc + 1 )||/ / ) 
x ^iV"' 0-1 (nilinjlj,ni'limjilf) + A r ' c_1 (njljnili,nj>ljirii>li>)J . 

(43) 

The rest of two-particle Breit-Pauli operators that we did not investigate 
so far are the two-body Darwin and spin-spin-contact terms. They do not 
bring any additional difficulties into the investigation of Hamiltonian, but 
for the sake of completeness of presentation we will discuss them briefly 

The two-body Darwin operator Hd 2 (see for more detail Nikitin and 
Rudzikas £ 4J, as well as the spin- spin- contact operator Hssc (see Shalit and 
Talmi |40j and Feneuille |41j). both have the following tensorial structure: 

H = y\tf(fcfco,ooo) (44) 

k 

These two terms are included into calculation by adding to the radial integral 
Rk (nih^'h' 1 n jlj n j'lj') a term 

(2k + 1) X (niiiUi'li/ , UjljUj/lj') , 

where 

X (nilinjlj,niiliinjilj/) 
= i fo°° Io°° P fa) Pj fa) %6 fa - r 2 ) P V fa) P r fa) d n dr 2 . (45) 

The expression (j!8|) has a series of terms, and thus at a first glance seems 
to be difficult to apply. For this purpose in the next sections we shall discuss 
in more detail: 
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• Obtaining the expression by Wick's theorem. 



• The compact written form of all terms, using the extended graphical 
technique. 

• Obtaining the values of recoupling matrix and of the standard quanti- 
ties. We shall also compare the existing methodologies of finding an- 
gular parts, showing the advantages and shortcomings of one method- 
ology or another. 

This methodology is fully applicable also for the one-particle operator 
(see Gaigalas and Rudzikas [22]) Gaigalas et al [HH]). As this operator does 
not cause big problems in the atomic physics, we will not stop for its details. 

3 Wick's theorem 

Wick's theorem in the second quantization formalism is formulated as follows 
(see Wick 02]; Bogoliubov and Shirkov [IS]): If A is a product of creation 
and annihilation operators, then 



where {^4} represents the normal form of A and |^4j represents the sum 
of the normal-ordered terms obtained by making all possible single, double, 
... contractions within A. Based on the identification in Bogoliubov and 
Shirkov |43| , the operator is presented in normal form when all of the opera- 
tors of annihilation included in it are transferred to the right of the creation 
operators. 

Usually, Wick's theorem is applied when treating complex operators that 
are represented by a large number of second quantization operators in a non- 
normal product form. In atomic physics such operators are used in perturba- 
tion theory (see Fetter and Walecka [H] ; Lindgren and Morrison [2^ ; Merke- 
lis et al jUJ) and in the orthogonal operator method (see Uylings Most 
often the Wick's theorem is applied to the products of second quantization 
operators that are not tensorially coupled (see Lindgren and Morrison |21j^. 
While applying the perturbation theory in an extended model space, two 
different groups of second quantization operators are defined (see for details 
in Lindgren and Morrison |2J). The second quantization operators acting 
upon core states belong to one group, whereas the operators acting upon 
open and excited shells belong to another one. These two groups are very 
different in applying Wick's theorem to them. 




(46) 
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In the first group, the creation operators are re-named to annihilation 
operators and are called the hole absorption operators, while the annihilation 
operators are re-named to creation operators and are called the hole creation 
operators. The creation operators of the second group are called the particle 
creation operators, while annihilation ones are called the particle absorption 
operators. Such a division of second quantization operators into two groups 
is called the particle-hole formalism. 

Merkelis et al [22] have proposed to use the so-called graphical ana- 
logue of Wick's theorem in perturbation theory (see Gaigalas et al |45j . 
Gaigalas [SI)- It is an efficient tool for obtaining the normal products of 
second quantization operators in a coupled tensorial form. Before applying 
this theorem, particular second quantization terms are in a normal product 
in coupled form. In addition, this theorem is applied in the particle-hole 
formalism, too (see Merkelis et al [22], Merkelis et al |47j ) . 

In all the cases mentioned above, the Wick's theorem is applied for the 
most general case of operators, i.e. the shells that are acted upon are not 
detailed. But, in the case of the extended model space, the group that 
the second quantization operators belong to, depending on the electronic 
structure of atom or ion under investigation, is defined. 

Wick's theorem is not applied in investigations of ordinary physical op- 
erators. Gaigalas et al [33] proposed a new version of the Wick's theorem 
application, where the optimal tensorial expression of any two-particle op- 
erator is easily obtained. The specificity of Wick's theorem application in 
this case lies in applying it only when the shells that are acted upon by 
the secondary quantization operators are known, i.e. it is applied for each 
term G(ij, i'j') separately. Such an interpretation of Wick's theorem bears 
similarity with the particle-hole formalism. The only difference is that in 
this case the second quantization operators are differentiated formally not 
on the basis of structure of atom under investigation, but on the basis of 
shells acted upon. 

This is done in the following way. The second quantization operators 
acting upon a shell with a lowest index are attributed to the first group. 
Those acting upon a shell with a next-lowest index are attributed to the 
second group, etc. In the most general case we have four distinct groups. 

For example, suppose we have a two-particle operator G(25,32), 

(5(25,32) = \ a 2 a 5 a\al (2,5|g|3,2), (47) 

(where 2 = ri2l2smi 2 m S2 , 3 = nsl3smi 3 m S3 , 5 = n^smi 5 m Sh ), the matrix 
element of which between the functions 



18 



^5 {LSM L M S ) = (n 1 hn2l2L 1 2Si2n 3 l 3 Li23Si23n4kLi234:Si23in5hLS\ and 
(L'S'M' L M' S ) = |nl/ 1 n 2 / 2 L / 12 5i 2 n3^3^i23«5^23«4/4^ / l 2 34^234™5/5^ , S , ) we 
must obtain. Then the operators acting upon the second shell are attributed 
to the first group, the ones acting upon the third shell - to the second, and 
the ones acting upon the fifth - to the third group. 

Assuming that all the operators from first group are the creation ones, 
and the rest are annihilation operators, we apply the Wick's theorem. Thus 
we obtain that the operators from the first group are one beside another, 
and all are positioned after the operators from the first group. Thus we 
obtain: 

— a2a\a^a\. (48) 

After that, we apply Wick's theorem assuming that the operators from 
the first and the second groups are creation ones, and the rest of them are 
annihilation operators. Thus we obtain that the operators from the second 
group are one beside another, and all are positioned after the first group 
operators: 

a2a\a\a^. (49) 

If in the product that we investigate there are operators of second quan- 
tization acting upon four distinct shells, then we apply Wick's theorem once 
again, assuming that operators from the first, second and third groups are 
creation ones, and from the fourth group - annihilation ones. In this case the 
Wick's theorem is applied to the second quantization operators in uncoupled 
form. 

Prom 1)11)1 we see that in second quantization a two-particle operator is 
written as a sum, where parameters i, j, i! , f run over all possible arrays 
of quantum numbers. So, the greater the number of open shells in bra and 
ket functions, the greater the number of terms G(ij, i'j') in the expression of 
two-particle operator. It is obvious that all these terms must be systematized 
in order to obtain in general case the most efficient tensorial expression of a 
two-particle operator, in the way described above. 



19 



Table 1: Distributions of shells, upon which the second quantization oper- 
ators are acting, that appear in the submatrix elements of any two-particle 
operator, when bra and ket functions have u open shells (see Gaigalas et 
al 33 ) 



No. 


di 


aj 


4 


4 


submatrix element 


1. 


a 


a 


a 


a 


(...^^•••IIGK^^-n^n^;.)!!...^^...) 


2. 


a 


P 


a 


P 


(„.n a e...v^-ii 


3. 


P 


a 


P 


a 


4. 


a 


P 


P 


a 


G^iknjljn'jl'^l'j) 
\\...nj^...npl 13 ...) 


5. 


P 


a 


a 


P 


6. 


a 


a 


P 


P 


(...n Q e...n^ V ^..||G||...n a e-^..V^...) 


7. 


P 


a 


a 


a 




8. 


a 


P 


a 


a 




9. 


p 


P 


P 


a 


G^iknjljn'jl^n'jl'j) 
\\...nX<*+\..n f3 lfi-\..) 


10. 


p 


P 


a 


P 


11. 


p 


7 


a 


7 




12. 


7 


P 


7 


a 


. .TIq/Iq a Tl^l Q ^ Tl-yl-y ^ . . . | | 


13. 
14. 


7 
P 


P 
7 


a 

7 


7 
a 


Ginilirijljri^rijl'j) 
\\...nX^n4"-\l^...) 


15. 
16. 


7 
7 


7 
7 


a 

P 


P 
a 


(...n a l% a nplp f> n 1 ly' 1 ... ||G(njZjnjZj7i^n^) 
\\---n a l^ a ~^^npl n^yl-y 1 ...) 


17. 
18. 


a 
P 


P 
a 


7 
7 


7 
7 


(..^J^ a n^ n 1 l^\..\\G(n i l i n j l j n' i iy j r j ) 
\\...n a l^ a ^nplp 13 Uyl-y 1 ...) 


19. 


a 


P 


7 


5 




20. 


P 


a 


7 


5 


(n a la a np Ip n 7 1^ 1 n s lf s \\ 


21. 
22. 


a 
P 


P 
a 


5 
5 


7 
7 


Ginilinjljn'^n'j l'j ) 
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Table 1 (continued) 



No. 


CLi 


aj 


4 


4 


submatrix element 


23. 
24. 
25. 
26. 


7 
<5 


5 
5 

7 
7 


a 
/3 
a 

/? 


Q! 
P 

a 


(n n l Na nala P ru,fi"'nsl? t 1 1 

\' v Ot v OL P p 1 tj 'y j f ,tJ O v fo | 

(^(nilinjljn'^n'jlj) 


27. 
28. 
29. 
30. 


a 

Ct 

7 
7 


7 

T 
/ 

a 
a 


/? 
/? 


S 


P 
5 


(nJ^nRl^nJ^nd? 1 1 1 

V Q G P p 'fywy ,(j O v q \ 1 

G^riilinjljn'^n'jlj) 


31. 
32. 
33. 
34. 


/3 
5 

P 
5 


5 

P 
5 

P 


a 
i 

7 
a 


7 
a, 
a 

7 


{n r ,l^ a nFtl^ 13 n^l^ 1 nxl^ 6 1 1 

\ li Q P P I I O A {I 

Ginilinjljn'^rijlj) 
WnJ^n^n^nsl^- 1 ) 


35. 
36. 
37. 
38. 


a 
5 
a 
5 


5 
a 
5 
a 


/? 
7 
7 


7 
/9 
/3 

7 


GiniliTijljn'^n'jl'j) 
\\n a l^-Wf + \^ +1 ^ s - 1 ) 


39. 
40. 
41. 
42. 


P 
7 
P 
7 


7 
P 
7 
P 


a 
a 


5 
a 
a 
5 


(rial^nplp 13 rijl^ n s lf 5 
GiniliTijljn'^n'jl'j) 

iKe+v?~V7 7 ~Vf +i ) 



In the work by Gaigalas e£ aZ |HH] it is chosen an optimal number of 
distributions, which is necessary to obtain the matrix elements of any two- 
particle operator, when the bra and ket functions consist of arbitrary number 
of shells. This is presented in Table 1. We point out that for distributions 2-5 
and 19-42 the shells' sequence numbers a, P, 7, 5 (in bra and ket functions 
of a submatrix element) satisfy the condition a < P < 7 < 5, while for 
distributions 6-18 no conditions upon a, P, 7, 5 are imposed (This permits 
to reduce the number of distributions). For distributions 19-42 this condition 
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is imposed only for obtaining simple analytical expressions for the recoupling 
matrices R Xj,X' i} X'j,A bra , A ket , r) . This will be discussed in more detail 
in the next section. 

So, in the way that is described earlier, the Wick's theorem is applied, 
assuming that the second quantization operators acting upon shells a, f3, 7 
and 5 belong to different groups. 

In the next section we discuss the way to obtain irreducible tensorial 
form of these distributions. In addition, the arguments will be given in 
evidence of superiority of the obtained tensorial expressions against other 
expressions known in the literature. 

The methodology presented in this section demonstrates the way to ob- 
tain optimal arrangement of the second quantization operators, for any two- 
particle operator. It can be applied without restrictions for obtaining the 
optimal tensorial form of two-particle terms of orthogonal operators and of 
perturbation theory operators, too. 

4 Graphical methods for two-particle operator 

The graphical technique of angular momentum is widely used in the atomic 
physics: see Yutsis et al [3S], Jucys and Bandzaitis [35], Brink and Satcher [5D| . 
El-Baz It is applied efficiently both in the coordinate representation 

(see for example Jucys and Bandzaitis 49 ), and in the second quantiza- 
tion formalism (Gaigalas et al |52| ) . The use of it allows one to obtain the 
analytical expressions for the recoupling matrices conveniently (see for ex- 
ample Kaniauskas and Rudzikas 53 ), to investigate the tensorial products 
of operators (see for example Jucys et al 1.54), and to seek for the matrix 
elements of operators (see for example Huang and Starace [SSI)- Gaigalas 
and Merkelis f2] have proposed a graphical way to obtain the values of ma- 
trix elements when the operator is a many-particle (one-, two-, three-, etc.) 
one and has irreducible tensorial form. For example, when a two-particle 
operator is considered, it has the form (|12|) or (|14|). The matrix elements 
in this methodology are expressed not only in the terms of coefficients of 
fractional parentage or reduced coefficients of fractional parentage, but in 
terms of standard quantities U k and V kl , too. Gaigalas et al 33 have pro- 
posed to calculate the matrix elements by using the tensorial expressions 
for such two-particle operator, that take full advantage of Racah algebra. 
In this case the tensorial form of operator depends on the shells that the 
operator acts upon (distributions 1-42 from Table 1). This is the difference 
of this methodology from other. It is most convenient to obtain the tensorial 
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expressions for 42 distributions graphically, using the generalized graphical 
technique by Gaigalas et al L 52 . In such not only the similarities 

between different distributions are easily seen, but also the compact graph- 
ical representation of the obtained expressions is possible. We will stop for 
details on this in the present section. 

A two-particle operator may be represented graphically by a Feynman- 
Goldstone diagram D\ from Figure ^ (Lindgren and Morrison [2J| ) ■ As it 
is shown in the paper Bolotin et al jSH], the Feynman-Goldstone diagrams 
are topologically equivalent to the angular momentum graphs. Due to that, 
an irreducible tensorial form for every Feynman-Goldstone diagram may be 
obtained (see Merkelis et al |1Z|)- The graph D2 is the angular momentum 
graph corresponding to the diagram D\. So the two-particle operator will 
be written down as follows: 



G(ij,i'f) = D 1 = -i y, Efa'^i^r^x 

mx.mx m x ,m x , p 
J i j 

x (niAinjAjH^^^lln^A;) D 2 afejafe] a^Ja^, (50) 

i i 

It must be noted that in expression (|50~|) the projection m\ i of am x , , as 
well as that of the momentum line Aj in graph D2 are the same. This is 
also to be said about the remaining operators of second quantization and 
the three open lines af graph D2 ■ 

As it has been mentiond in section 2, this operator has two tensorial 

forms, (|12j) and These may be represented graphically, since the cre- 

r\ \ ~( A ') 
ation operator a ( l >, as well as operator a J respectivelly, are graphically 

denoted by diagrams -D3 and D4. 

The first form (|12[) of two-particle operator Gi(ij,i' f) is represented as: 



G I (ij,i'j') = D 1 = -l E(" 1 )^ P [«i.«2,*i> ^ 2 ] 1/2 x 



K12K 12 (T12(T 12 



la In 



.12 J 1Z 

/j lj K\2 



s 


s 


„l 






k 


s 


s 


0\2 



whereas the second (|13|) : 



(51) 
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riiXi 



KK. 



kk 



rijXj 



n>X> 



Di 



+ 
A' 



— I — 

kk 

D 2 



A' 



riiXi 



D 3 



rijXj 



n'X' 



kk 



rijXj 



niX t 
+ 



rijXj 



«12 &12 k' 12 o\o 

kk 



n'X 



n>X> 



D 5 



Da 




D 7 D 8 D 9 



Figure 1: Diagrams for two-particle operators. 
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Gii(tj,i'f) = D 5 + D 7 = 

= ^E(- 1 ) fe ~ P (n^njX^^^lKX^X'j) x 
p 

^ ? * }<5(n^,n^)^9}. (52) 

We emphasize here that the winding line of interaction in the Feynman- 
Goldstone diagram corresponds to the operators of second quantization in 
the normal order (Figure ^ D\). Whereas the dotted interaction line indi- 
cates that the second quantization operators are ordered in pairs of creation- 
annihilation operators. In the latter case first comes the pair on the left side 
of a Feynman-Goldstone diagram (Figure D§). Such a notation of two 
kinds for an interaction line is meaningful only for two-particle (or more) 
operators, since for any one-particle operator both the winding and dotted 
lines correspond to the same order of creation and annihilation operators. 

From expressions ([51 jl . (|52|) we see that the two-particle operator in the 
first form is represented by one Feynman-Goldstone diagram Di, whereas 
in the second - by two diagrams D§ and Dj, The diagrams, corresponding 
to tensorial product, have the following algebraic expressions: 




Mi 



x a y 



(K12CT12) 



x a 



x a 



( A 



(ftlCTl) 



a (Aj) x a 



(K2CT2) 



(kk) 

P-P 
(kk) 

P-P 



Da 



a (Al) x a v ' J 



(kk) 



P-P 



(53) 



(54) 



(55) 



The positions of the second quantization operators in the diagrams Dq, 
and Dg define their order in tensorial products: the first place in ten- 
sorial product occupies the upper right second quantization operator, the 
second - lower right, after them the upper left and lower left operators follow. 
The angular momenta diagram defines their coupling scheme into tensorial 
product. For more detail see Gaigalas and Rudzikas [32| . 
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Figure 2: Diagrams for distribution 7, f3, a, 7. 



As it has been mentioned earlier, these two forms do not always take 
full advantage of the Racah algebra (see Gaigalas and Rudzikas 32 ) . The 
expression (|18|) has no such shortcomings. Now we will demonstrate the 
way to obtain graphically a tensorial expressions for particular distributions 
1-42 from Table 1. 

We take the distribution 7, /3, a, 7 (13 form Table 1) as an example for 
investigation. Then the Feynman-Goldstone diagram of operator G(7/3, cry) 
is -D10, the angular momentum graph is Du, and the second quantization 
operators are in the following order: 

a^apa^a^. (56) 
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Applying the Wick's theorem as described in Section 3, and assuming that 
the operators acting upon shell 7 belong to the first group, the once acting 
upon (3 belong to the second, and the ones acting upon a belong to the third 
group, we obtain the following order of operators: 

a^apa^. (57) 

Using the generalized graphical technique of angular momentum by Gaigalas 
et al we couple the operators of second quantization into tensorial 
product D12 (see Figure |3J): 



D 



12 



(K12012) 



4^ 



x a 



.(A a ) 



(kk) 



P~P 



(58) 



In the course of obtaining D\2 graphically, a recoupling matrix D13 appears, 
whose analytical expression is readily obtained by using the graphical tech- 
nique of Jucys and Bandzaitis . All the needed expressions are obtained 
in the same way. 

In obtaining these expressions, as well as in representing them graphi- 
cally, it is very convenient to use the rule of changing the sign of a node, 
existing in the graphical technique of angular momentum (see Jucys and 
Bandzaitis [HI]). Therefore now we will treat an example of using such a 
rule. 

Suppose, we have the following correspondence between diagrams (Fig- 
ure EJ): 



D 14 — ► D 1S , (59) 
in which the second quantization operators are in the order 

a (A 3 ) S (A 4 ) q (Ai) 

a( X2 \ Our goal is to obtain the diagram corresponding to the order a^ Al ^ 
a(A 2 ) a (A 3 ) S (A 4 )_ 

Bearing in mind that the second quantization operators 
anticommute with each other and they all act on different electronic shells 
and we are not changing the order of their coupling into tensorial product, 
we arrive at 

D U — > (-l) 4 £ie = D w . (60) 

Let us also discuss another situation: we have defined the disposition 
of the operators and we want to couple them into certain tensorial prod- 
uct. Suppose that we want to represent graphically the following tensorial 
product: 
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Figure 3: Diagrams for graphical transformations. 
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(ftlCTl) 




(61) 



a 



x a 



For this purpose we have to rearrange the generalized Clebsch-Gordan 
coefficient, which is defining the scheme of coupling of the operators into the 
tensorial product. It is easy to notice that this coefficient will properly define 
the tensorial product, if we change the sign of the vertex "a" in diagram 



The procedures described are fairly simple, however, they are sufficient 
for the majority of cases. The more complete description of this general- 
ized graphical approach may be found in Gaigalas et al (13]) Gaigalas [15] . 
Gaigalas and Merkelis 52 . 

All the analytical expressions for distributions 1-42 from Table 1 are 
presented in the paper Gaigalas et al [22| . They are written down using the 
generalized graphical methodology of angular momentum, and the vortex 
sign change rule, which was discussed in this section. As a consequence 
of that, the analytical expressions for 42 terms may be written down via 6 
different expressions. This, undoubtedly, facilitates a lot the implementation 
of methodology proposed in Gaigalas et al 33 . 



5 Matrix Elements Between Complex Configura- 



In this section we will discuss several ways to obtain matrix elements of 
a two-particle operator. As it was mentioned earlier, up to now the Fano 
calculation scheme ^U] is the most popular one. Its general expression when 
a two-particle operator acts upon different shells is presented in 0. 

The general expression for a matrix element in other cases is similar. For 
example, when the operator acts only upon one shell, we have 




(62) 



tions 
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(LS) \\G^ K2k ^ a '^\\^ et (US')) ~ 
Ni(Ni-l)x 

{T} 

xi? d ^Ai, A;, Aj, Ai,Ki,K 2 ,K,ai,CT2,cr, A fera , A fcei ) x 
x (niAiniAill^^'^^HniAiniAi) . 

(63) 

As we see here, in contrast to (JSJ, the summation is performed only over 
one array n, A of quantum numbers, because the operator acts only upon one 
shell. But here a summation over {T} occurs, though, which indicates the 
summation over arrays of intermediate terms a"L'-S", a'-'L'-'S'/', a'^ L'^ S'^ . 

Remembering the relationship between a coefficient of fractional parent- 
age and a reduced matrix element of a second quantization operator (see 
Spakauskas et al (57], Rudzikas and Kaniauskas 20 ): 

(l N aLS\\aW\\l N - 1 a'L'S') = 
= (-i)N+(N+lMN) ^n[l^] (l N aLS\\l N ~ l {a'L'S 1 ) , l) (64) 

we see that the Racah algebra in expressions Q, (|63(l is used only on the 
level of coefficients of fractional parentage. In separate cases, e.g. when 
the two-particle operator acts upon one or two shells, it is possible to use 
expressions which exploit the Racah algebra at a higher level, i.e. to take 
more advantage of the tensor algebra elements (see Judd |27J, Jucys and 
Savukynas 3 ). For example, let us investigate the case when a matrix 
element is calculated for bra and ket functions having one shell only. The 
tensorial forms l|12[) and (|14|) are of value here. Taking the second one of 
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these, we have 



(nl N aLS\G n \nl N a'L'S') = 
\ (nXnX\\g^ lK2k ' aia2 ^\\nXnX^ {[ki,k 2 ,<ti,o >\ 



x{nl N aLS\ 



x a 
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(kicti) 
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x(nl N aLS\ 
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-1/2 



\nl N a'L'S') 



\nl N aL'S') 



(65) 

Using the relationships between the tensorial product of creation and an- 
nihilation operators and the tensorial quantities U k and V kl (see Rudzikas 
and Kaniauskas [2111); the expression (|65|) for matrix elements can be writen 
down in terms of U k and V kl . In comparing (|63j) to (|65[) we see that the 
summation over intermediate terms a'{L'{S'/, a'fL'fS'/', ol^ L'^ S'^ is al- 
ready performed in expression ()65|) . So, in this case the Racah algebra 
is exploited at the level of standard quantities U k ir V kl . This simplifies 
calculations a lot: 

• For zero matrix elements are easily tracked down from triangular con- 
ditions even before the actual calculation is performed. In case 
only the triangular conditions 5(L,k,L') and 5(S,k,S') are present, 
but their number may be greater in other cases. (In the above, the no- 
tation 5(L, k, L') means the triangular condition | L—L' \ < k < L+V .) 

• The tables of standard quantities (see Nielson and Koster _SB; , Karazija 
et al^, Cowan 0) may be used. 

• The recoupling matrix is simpler in this case, and it has an analytical 
expression. 

So, the expressions exploiting the Racah algebra at the level of U k and 
V kl are much more advantageous than Q. Such expressions are obtained 
for all physical operators. For example, the expressions for spin-other-orbit 
operator are presented in papers Horie [M|> Karazija et al (HJ and Vizbaraite 
et al the ones for spin-spin operator - in papers Horie 60^ and Karazija 
et al and the ones for orbit-orbit operator in the monograph Jucys and 
Savukynas 3 . The shortcoming of the expressions of this type is that the 
Racah algebra is exploited to its full extent in separate cases only. This is 
discussed in detail in paper by Gaigalas et al [3*2] . 
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Gaigalas et al [33] have proposed a methodology which allows one to take 
all the advantages of the Racah algebra in the most general case. According 
to the approach by Gaigalas et al [33], a general expression of submatrix 
element for any two-particle operator between functions with u open shells 
can be written down as follows: 

(V^ ra {LS) \\G^ k ' a ^\\^ ket {US')) = 
E ^u a (LS)\\G (n^n^Mitn'ti) (L'S')) = 

E E Y,(- 1 ) A@ '{^ n ^v<K,n' 1 X' J ,E)x 

xT (n i X i ,n j X j ,n' i X' l ,n' j X' 1 ,A bra , A ket , 3, r) x 
x J R(A i ,A J ,A^,A;.,A b ™,A fce4 ,r) . 

(66) 

In calculating the spin-angular part of a submatrix element using (|66|) . 
one has to compute the following quantities (for more detail see Gaigalas |33|): 

1. The recoup ling matrix R ( Xi , Xj , X'i , X'j , A bra , A ket , r) . This recoupling 
matrix accounts for the change in going from matrix element 

{xp b u ra (LS) \\G\\ip ket {L'S')), which has u open shells in the bra and ket 
functions, to the submatrix element 

T (n i Xi,n j X j ,n' i X' i ,n' j X' j ,A bra ,A ket , S, r) , which has only the shells be- 
ing acted upon by the two-particle operator in its bra and ket func- 
tions. 

2. The submatrix element T (riiXi, rij Xj , n[ A • , n'j Xj , A hra , A ket , E , Tj , which 
denotes the submatrix elements of operators - (j23|) . Here T refers 
to the array of coupling parameters connecting the recoupling matrix 
R (Ai, Xj, X' i ,X'j,A bra , A ket , r) to the submatrix element. 

3. Phase factor A (for more detail see Gaigalas 33 ). 

4. 0' ^niXi,njXj,n' i X' i ,n'jX'j,'B \ , which is proportional to the radial part 
and corresponds to one of G (nA, E),...,0 (n Q X a , npXp, n 7 A 7 , n$Xs, S). 
It consists of a submatrix element (niXiUjXj \\g( K ^ K i k ' a ^ k ) | ^[X'jTi'jX'A , 
and in some cases of simple factors and 3nj-coemcients (for more detail 
see Gaigalas [33]). 
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In the next section we shall discuss finding of the recoupling matrix. Now 
we shall analyse the submatrix element T frij Aj, nj\j,n' i X' i ,n'jX'j,A bra ,A ket , H, I 
As the angular part of expression contains the tensors ifTU]) - 1)23(1 . so 
we will discuss the derivation of submatrix elements of these operators, and 
present the expressions for these quantities. It is worth noting that these 
tensorial quantities all act upon the same shell. So, all the advantages of 
tensor algebra and the quasispin formalism may be exploited efficiently. 

We obtain the submatrix elements of operator (|19|) by straightforwardly 
using the Wigner-Eckart theorem in quasispin space (see Rudzikas 0): 



(l N aQLS\\a$ q s) \\l N ' a'Q'L'S' 



[Q] 



-1/2 



Q' 1/2 Q 



M' r 



x (l aQLS\\\a^\\\l a'Q'L'S' 



m q Mq 



(67) 



where the last multiplier in (ffiTj) is the so-called completely reduced (reduced 
in the quasispin, orbital and spin spaces) matrix element. 

The value of the submatrix element of operator 1120(1 is obtained by 



nl N aQLS\\ 

E wr 1/a 



a {qX) x a {gX) 



(feifea) 



nl N ' a'Q'L'S 



q q e 

m q i m q2 m € 



Q' 



Q 



Mq m e M Q 



x fnl aQLS\\\W^ k ^\\\nl a'Q'L'S' 



(68) 



On the right-hand side of equations (|o7|) and (|B^)l only the Clebsch- 
Q' e Q 



Gordan coefficient 



M Q 



m f Mr 



depends on the number N of equiva- 



lent electrons. 

(nl aQLS\\\W^ k ^\\\nl a'Q'L'S'^j denotes reduced in quasispin space 
submatrix element (completely reduced matrix element) of the triple ten- 

(ek 1 k 2 ) 



sor 



(qls) x a (qls) 



It is related to the completely 



reduced coefficients (subcoefficients) of fractional parentage in a following 
way: 
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X 



(nl aQLS\\\W^ klk ^\\\nl a'Q'L' S' 
= (_ 1 jQ+L+S+Q>+L>+S'+E+k 1 +k 2 ^ ^ k 2 ] 1/2 

x (laQLS\\\a {qls) \\\la"Q"L"S")x 

a"Q"L"S" 

x (l a ,/ Q' / ^ // fi ,,/ |||o ( «^|||i a'Q'L' S') x 
q q e \ J I I k\ 1 J s s k 2 
Q' Q Q" \\ L' L L" \\S'S S" 



(69) 



In the other three cases (|21|). (|22|) . (|23j) we obtain the submatrix elements 
of these operators by using (2.28) of Jucys and Savukynas [3]: 

(n/^ aQLS|| [f( KiCTi ) (nA) x G( K2<T2 ) (nA)] m \\nl N ' a'Q'L' S') = 

= i _ 1) L+S + L>+S'+2k [k]x 

x J2 ( nlN aQLS\\F^ K1 ^ (nA) | \nl N " a"Q"L"S") x 

a"Q"L"S" 

x(nl N " a"Q"L"S"\\G^^ (nX)\\nl N ' a'Q'L' S')x 

Kl K2 k 1 I 0"! 0"2 1 



L' L L" f|S" S 5" 



(70) 



where (nA), G( K2ff2 ) (nA) are one of or (1201) and the submatrix 

elements correspondingly are defined by (|67jl. (|68|) and H69|) , A^" is defined by 
the second quantization operators occurring in 

( n A) and G^* 7 ^ (nA). 

As is seen, by using this approach Gaigalas et al j^Hl, the calculation of 
the angular parts of matrix elements between functions with u open shells 
is reduced to requiring the reduced coefficients of fractional parentage or 
the tensors (for example W^ eklk2 ^ (nl,nl)), which are independent of the 
occupation number of the shell and are acting on one shell of equivalent 
electrons. 

The main advantage of this approach is that the standard data tables 
in such a case will be much smaller in comparison with tables of the usual 
coefficients U k , V klk2 (see Jucys and Savukynas [3]) and, therefore, many 
summations will be less time-consuming. Also one can see that in such an 
approach the submatrix elements of standard tensors and subcoefficients 
of fractional parentage actually can be treated in a uniform way as they 
all are the completely reduced matrix elements of the second quantization 
operators. Hence, all methodology of calculation of matrix elements will be 
much more universal in comparison with the traditional one (see Cowan [2], 
Jucys and Savukynas 0, Wybourne [HI]). 
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6 Recoupling Matrix 



While seeking the matrix elements of one- or two- particle operators, it is 
necessary to obtain the values of a recoupling matrix 

we use the methodology by Gaigalas et 
al [33] ( see expresion (|66|l). or the recoupling matrices 

R d ^Aj, Xj, \l, \p K%, K2, k,gi,<J2,<j, A bra ,A ket ^j and 

Xj, X[, Xj, k\, K2, k, ax, (J2, cr, A bra , A ket ^ ; if we use the methodol- 
ogy by Fano ^DJ (see expresion @). In the case of several open shells the 
expressions for matrix elements of every physical operator are published, 
where the recoupling matrices are in the form of simple factors. Usually all 
these are written in the coordinate representation. They can be found in 
Karazija et al |2DJ, Jucys and Savukynas [3J, Karazija jDS], Rudzikas [2D] . 

Meanwhile, for the more complex configurations, i.e. the ones having 
many shells, the recoupling matrices are much more complicated. Beside 
that, the complexity of a two-particle operator adds to this. When the ten- 
sorial structure of an operator is complex, the recoupling matrix is rather 
complex, too, e.g. the spin-other-orbit operator (see Gaigalas (M]). While 
attempting to calculate the angular part of matrix elements in all the men- 
tioned general methodology for calculating the recoupling matrices 
is necessary. It has to be efficient, too, because the speed of calculation of 
angular parts of matrix elements depends on that. 

The majority of methodologies to obtain angular parts are based on 
the Fano ^Dj calculation scheme (see for example Glass |DJ, Glass and Hi- 
bbert |7j, Grant [S]). In finding the matrix elements using this, one of the 
tasks is to obtain the recoupling matrices for direct and exchange terms. 

Let us treat a matrix element of a two-particle operator acting upon four 
distinct shells 

(LSM L M S ) \\G(ij,i'j')\\^ et (L'S'M' L M' S )) . (71) 

For such an operator, the recoupling matrix in L - space, using graphical 
technique of Jucys and Bandzaitis [ID], is represented in Figure [1J As we 
speak of calculations in LS - coupling, the analogous recoupling matrix in 
S - space has to be calculated. In the Figure the tree containing B nodes 
represents the bra function, and the one with K nodes represents the ket 
function. Bra, as well as ket function, both contain u open shells in this 
case. 

Now let us treat a case when three distinct shells are acted upon by the 
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Kj K s Kj K d Kj K„ D 



Figure 4: The recoupling matrix for the Fano angular calculation scheme. 

operator. For example, the matrix element of operator G(ij,ij') is sought. 
Then in the recoupling matrix the node III will be on Lj line (see Figure 0}. 
This means that there will be two nodes on the resulting orbital moment Li 
line of shell i. These will be nodes I and III. The resulting orbital moment 
Li of the shell i of bra function is attributed to the line connecting nodes 
Bi and I. The resulting orbital moment L\ of the shell i of ket function is 
attributed to the line between nodes III and K^. Finally, characteristics 
L'l of all possible intermediate momenta for that particular shell, which are 
summed over, are attributed to the line between nodes I and III. So, in 
treating the operator acting upon three distinct shells, an intermediate sum 
in the recoupling matrix appears. Similarly, two intermediate sums occur 
when two distinct shells are acted upon by the operator. And when it acts 
upon one shell only, three intermediate sums are present. The lines k, k' and 
K in the recoupling matrix represent the structure of two- particle operator 
in / - space. See Tutlys (SHI, Grant [S] for details. 

This is also valid for the exchange term, only the line 1^ must be con- 
nected to the node D, and line ly - to the node C. 

The first program to calculate the recoupling matrices of this type, 
NJSYM, was written by Burke ^]. It performs the calculations in two 
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stages: 1) the recoupling matrix is expressed as a sum of products of the 6j 

- coefficients; 2) this expression is used in calculation. 

Tutlys wrote a program to calculate angular parts of matrix ele- 
ments, ANGULA, which expressed the recoupling matrix in terms of Clebsch 

- Gordan coefficients before the actual calculations. While finding the recou- 
pling matrix by the Clebsch - Gordan coefficient summation, this program 
eliminates trivial coefficients from the expression. 

Bar-Shalom and Klapisch [^j developed a new program NJGRAF. This 
program calculates the recoupling matrix in several stages. On the basis of 
graphical methodology by Yutsis, Levinson and Vanagas j3H]) the recoupling 
matrix is analysed graphically and an optimal expression is found. After- 
wards, the value of recoupling matrix itself is calculated. An analogous 
program RECOUP was written by Lima [HI], and a program NEWGRAPH 
was written by Fack et al |BBj. All these (NJGRAF, RECOUP and NEW- 
GRAPH) are based on the same principle. An optimal analytical expression 
for the recoupling matrix is obtained by graphical method, and then the 
calculations are carried out according to it. But the optimal expressions 
they find are different quite often, and are not really optimal. 

As it was mentioned above, the methodology of angular calculation based 
on the Fano calculation scheme has a shortcoming that the intermediate 
sums appear in complex recoupling matrices. Due to these summations 
and the complexity of the recoupling matrix itself, the associated computer 
codes become rather time consuming. A solution to this problem was found 
by Burke et al They tabulated separate standard parts of recoupling 
matrices along with coefficients of fractional parentage at the beginning of 
a calculation and then used them later to calculate the coefficients needed. 

Computer codes by Glass Glass and Hibbert [7|, Burke et al [S], 
Fischer |69| . Fischer [7(J| and Dyall et al |71j utilize the program NJSYM 
(Burke JB]) or NJGRAF (Bar-Shalom and Klapisch for the calculation 
of recoupling matrices. Both are rather time consuming when calculating 
matrix elements of complex operators or electronic configurations when cal- 
culating matrix elements of complex operators or electronic configurations 
with many open shells. In order to simplify the calculations, Cowan [2] sug- 
gested grouping matrix elements into 'classes' (see Cowan Figure 13-15). 
Unfortunately this approach was not generalized to all two-electron opera- 
tors. Perhaps this is the reason why Cowan's approach is not widely used 
although the program itself, based on this approach is widely used. 

Gaigalas et al |33| proposed a methodology where the analytical expres- 
sions for recoupling matrices are obtained for the most general case. In 
this methodology, analogically as in Cowan [2j, the matrix elements are at- 
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Figure 5: Diagram Ri representing the recoupling matrix 
fl^A^A^A^A 6 ™, A* e * , T^j when a two-particle operator acts upon 
a single shell. 

tributed to four different groups. The operators acting upon only one shell 
belong to the first group (distribution 1 from Table 1), the ones acting upon 
two - to the second (distributions 2-10 from Table 1), upon three - to the 
third (distributions 11 - 18 from Table 1), and upon four - to the fourth 
group (distributions 19 - 42 from Table 1) respectively. Each group has a 
different recoupling matrix, R\ - i?2- They all are shown in Figures El - |H1 

6.1 One interacting shell 

Let us treat the recoupling matrix R\ first (see Figure EJ) ■ It represents the 
case when a physical two-particle operator acts upon a single shell. Here, 
similarly as in Figure^ the graphical technique of Jucys and Bandzaitis 
is used. The bra and ket functions are represented as in Figure EJ As the 
two-particle operator acts upon only one shell in this case, this operator 
may be treated as a single - particle one from the recoupling matrix point of 
view. Using a graphical rule allowing one to cut two lines and connect the 
loose ends, we disconnect the nodes B\ and K\ from the general recoupling 
matrix. The part separated from recoupling corresponds to delta - functions 
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Figure 6: Diagram R2 representing the recoupling matrix 
R(Xi,\j, K , X'j , A bra , A ket , r) when a two-particle operator acts upon 
a two shell. 
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Figure 7: Diagram i?3 representing the recoupling matrix 
R(Xi,\j, K , X'j , A bra , A ket , r) when a two-particle operator acts upon 
a three shell. 
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Figure 8: Diagram i?4 representing the recoupling matrix 
R(Xi,\j, K , X'j , A bra , A ket , r) when a two-particle operator acts upon 
a four shells. 
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I J 2 



J J2 4 



Figure 9: Diagram corresponding to the notation A\(+ — 
H — ; J1J2J12; J^JJ^)- 



only, 5(L\2, L' 12 ), S(L±, L2, L12). In addition, using graphical technique it is 
possible to cut all the nodes until Bi and K[ out of the general recoupling 
matrix. All the cut nodes contribute the same delta functions which may 
be written as <5(L 12 .. i __ 1 , L' 12 5(L 12 ,a-i, L h L 1%A+1 ). 

In the next stage, it remains to treat the part of recoupling matrix from 
the nodes Bi and Fi up to B u and K u . Beside that, the lines Li2..i-i 
and L' 12 j_i are connected in this remaining diagram. Before obtaining 
the analytical expression for the remaining recoupling matrix we introduce 

a notation A\{-\ 1 — ; J1J2J12', •HJJ^)- It will help us to describe the 

procedure of obtaining analytical expression for the remaining recoupling 
matrix in a simpler way. So, the diagram in Figure |§] is denoted as A\{-\ — 
H — ; J1J2J12; J^JJw)- It equals to 



^i(+-+-;JiJ 2 Ji 2 ;J3JJ32) = (-i) 2 - ? '[ji2,j32] 1/2 | J j\ J J ^ j 



(72) 



Now, using the graphical rule of cutting three lines, we cut out the 
nodes /, Bi and Ki from the recoupling matrix. Thus we obtain a diagram 

A\{ h+; KLjLi] Li2..i-\LL' 12 ..%)• This coefficient in the paper by Gaigalas 

et al is denoted as C±, and equals to 



c 1 = (-ir[L a ,rf/ 2 | k j L * ^ j 



(73) 



The diagram is cut from the nodes Bj and Kj up to B u and K u in the 
same way. It is expressed in terms of diagrams 
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Figure 10: Diagram corresponding to the notation A2(— + + — 
; J1J2J12', J3J4J34', Ji3,J2aJ)- 



^raax 

C2 {K, k m i n , k max ) = Yl 1 — ; KL' 12 A Li2...i; LiLi2...i-iL' 12 ...i-i)- 

(74) 

It is easy to notice that C2 (K, k m \ a , /c max ) corresponds to the coefficient 
C2 (K, fcmin , k max ) defined in Gaigalas et al [33] . 

6.2 Two interacting shells 

Before going into investigation of the recoupling matrix pictured in Figure 
we denote the diagram A2 (see Figure E3 ) as 

A2 ( — h H ; J1J2J12', «/i3«/24«/)- Its analytical expression is 

!Ji J2 J12 
J3 Ja J34 
J\3 J24 J 

(75) 

In the A2 diagram, first the signs of nodes 1-6, and then the moments are 
presented. 
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While investigating the recoupling matrix (see Figure |SJ) it is easy to 
notice that there is an additional graphic element. After cutting this out, 
we get a diagram the further investigation of which is very similar to the 
investigation of R\. That additional matrix element is 

A 2 (—+-\ ; L' X2 j_ x KxL V2 ,..j-\\ J'jK 2 Lj-, L' 12 -KL 12 ..j)- It is obtained 

by cutting the nodes Bj, II, III, Kj out of the diagram R 2 . After doing this, 
the diagram under investigation is split into two diagrams closely resembling 
the recoupling matrix R\. As they are cut in the very same way as R\ is, 
we will not stop at recoupling matrix R 2 for details. 

6.3 Three interacting shells 

In the methodology presented in Gaigalas et al 33], the recoupling matrix is 
represented by a diagram i?3 (see Figure EJ) when the two-particle operator 
is acting upon three shells. Similarly as for the recoupling matrix R 2 , the 
cutting out of nodes Bj, II, III and Kj leads to a diagram A 2 { — h H — 

; L' l2 j-\R~\B\ 2 ..j-\; J'jK 2 Lj; L' 12 jK?,L\ 2 ..j), and the cutting out of nodes 

B[, IV, V ir K[ leads to a diagram 

M{ — H ; L' 12 i ,_ 1 i^ 3 L 1 2..j'_i; J[,K^Lii\ L' 12 ^,KLv2..i>)- In that case 

the diagram splits up into three simple diagrams, the further investigation 
of which is analogous to that of R\ recoupling matrix. 

6.4 Four interacting shells 

The most complex recoupling matrix R^ (see Figure |HJ) comes into view 
when all the second quantization operators act upon different shells. The 
part of this recoupling matrix from nodes B\ and K\ to nodes B^i and K\i , 
and also from nodes Bji and Kj/ to nodes B u and K u is very similar to 
the recoupling matrices described above, therefore their investigation is the 
same as that of the diagrams described earlier. We will investigate in more 
detail the part of diagram from nodes B^ and Ky to nodes Bji and Kj/. 
That is pictured in Figure 1111 After cutting the diagram across four lines 
^12.. K^, K% and L' 12 v and connecting these according the graphical 
technique of Jucys and Bandzaitis [IHj with a generalized Clebsch-Gordan 
coefficient, then connecting the lines Li 2 ...i'-i, K3 and L' 12 with an 
ordinary Clebsch-Gordan coefficient, we obtain a diagram which is expressed 
via two 6j-coefficients. The exact analytical expression of this diagram is 
presented in the paper by Gaigalas et al (see expresion (29)). 

After cutting the nodes Bji, V, VI, VII and Kj/ out of the remain- 
ing part and investigating the obtained diagram according to the Jucys 
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Figure 11: The most complex part of the recoupling matrix R^. 



and Bandzaitis technique, we obtain a diagram which corresponds to a 
12j- coefficient. The analytical expression of that is presented in paper by 
Gaigalas et al [33], expr. (31). 

If a shell TLj/lji is beside shell n^Zj', i.e. there are no additional nodes 
between nodes Bp, K# and By, Kj>, then the additional graphical elements 
vanish. So, we obtained the full analytical expression of recoupling matrix 
-R4. If between the abovementioned shells a single additional shell is present, 
then we obtain an additional diagram which is expressed as a product of two 
6j-coefncients (see (30) from Gaigalas et al |33j ) . If the number of shells is 
greater, then the obatined diagram in the most general case is written via 
(32) from Gaigalas et al [3*H] . 



7 Conclusion 

The methodology that is based on the second quantization in coupled ten- 
sorial form, on the angular momentum theory in three spaces (orbital, spin 
and quasispin), on the Wick's theorem and on the generalized graphical 
technique of angular momentum, give the possibility to efficiently calcu- 
late the matrix elements of energy operators in general case (Gaigalas and 
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Rudzikas [52] . Gaigalas et al (33], this paper). The main ideas of this 
methodology are: 

1. A number of theoretical methods is known in atomic physics that fa- 
cilitate a lot the treatment of angular part of matrix elements. These 
are the theory of angular momentum, its graphical representation, 
the quasispin, the second quantization and its coupled tensorial form. 
But while treating the matrix elements of physical operators in gen- 
eral, the methods mentiond above are applied only some part, or very 
inefficiently 

An idea is presented and carried out in the work Gaigalas and Rudzikas |32| , 
Gaigalas et al |33| and this paper, of unifying all these methods in order 
to optimize the way to treat matrix elements of physical operators, to 
investigate efficiently even the most complex cases of atoms and ions. 

2. The Wick's theorem is widely known in the theory of atom. Up to now 
it was applied to the most general products of the operators of sec- 
ond quantization, i.e. when the particular arrays of quantum numbers 
for each operator of second quantization were not yet defined. In the 
methodology described there is proposed to apply the Wick's theorem 
for the products of second quantization operators where they have the 
values of quantum numbers already defined. This allows one, using 
the methodology of second quantization in coupled tensorial form, to 
abtain immediately the optimal tensorial expressions for any opera- 
tor. Then, in treting the matrix elements of physical operators, the 
advantages of a new modification of the Racah algebra are exploited 
to their full extent. 

3. While analysing the physical operators presented as products of ten- 
sors a^ qsl \ w( kqklka \ [a( qsl ^ xW^ kqklks ^ KqKlKa \ \W^ kqklk ^ xa^ qsl ^ KqKlKs \ 
^(kqhks) x wik'qklk'^KqKiKs) ^ - t j g pQggibig t obtain convenient an- 
alytical expressions for recoupling matrices that must be taken into 
account in finding the matrix elements (non-diagonal with respect to 
configuration included) of any physical operator between complex con- 
figurations (with any number of open shells). 

4. An idea is proposed and carried out, concerning the most efficient way 
to apply the tensor algebra and the quasispin formalism in the most 
general case, for the diagonal and the non-diagonal matrix elements 
as well, when the bra- and ket- functions have any number of open 
shells. 
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The combination of all these improvements allows to efficiently account 
for correlation efects practically for any atom or ion of periodical table. 
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